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Abstract 

Diffractive scattering of "fast" or "high energy " neutrons can give 
low energy nuclear recoils in the signal region for dark matter searches. 
We present a discussion using the 'black disc' model. This permits a 
simple and general, although approximate, description of this possible 
background. We note a number of its features. In particular there are 
mass number A dependent aspects which can be studied in setups where 
events on different nuclei are observable at the same time. These include 
the recoil energy distributions, and the A behavior of the cross section. 
We define a parameter which characterizes the recoil energy to be 
expected due to fast neutrons. It ranges from 100 keV on light nuclei to 
a few keV on heavy nuclei, and a general treatment is possible in terms 
of it, within the 'black disc' approximation. In addition the presence of 
inelastic processes on the nuclei at about the same level of elastic processes 
would be characteristic of fast neutrons. 

1 Introduction 

In the search for dark matter by direct detection in the laboratory, spectacular 
levels of background reduction have been achieved. In the "Commissioning Run" 
of the CRESST II detector, for example, the rate of posssibly interesting events 
was reported as 0.063 per kg-day [TJ, based on three events. For the about 600 
gm of sensitive material used in the dark matter analysis, this corresponds to less 
than one event per week. Similarly, the CDMS collaboration, after various cuts, 
had no events for 121 kg- days of germanium data .2. EDELWEISS discusses 
four events on the basis of 322 kg-days [3]. And most recently XENON100, with 
48 kg for 100 days, reports only three events [4]. 

With such low event rates even rare sources of background must be consid- 
ered. In particular the most worrisome background in underground experiments 
searching for dark matter through nuclear recoils are neutron-induced nuclear 
recoils. Neutrons scatter off nuclei and the recoils could closely resemble the 
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sought-for signal. Furthermore as we shall explain below, fast neutrons, which 
are difficult to shield, will give elastic nuclear recoils with the same general en- 
ergies as expected in the dark matter searches. Although estimates of neutrons 
from standard sources such as incoming muons or radioactivity in the surround- 
ing rock ( see the above references and for example jS]) give neutron fluxes which 
are quite small, the very good level of background suppression being reached 
implies that even very rare processes or perhaps unexpected contaminations 
must be considered. While there are elaborate Monte-Carlo program sets for 
carrying out standard background calculations, it must be said that they are 
not very transparent, and do not offer much insight into the nature of the effects 
under consideration. Hence it may be useful to have a discussion based on a 
simple physical picture, even if this is only semi-quantitatively accurate. 

Estimates generally lead to neutron sprectra dominated by low energies. 
Such fluxes can be strongly reduced by neutron shields. However such shields 
are typically hydrocarbon materials, which while very effective at low energy, 
become rapidly ineffective at higher energies. Thus the PDG tables [5] give the 
rather large value of 88 cm for the nuclear interaction length in polyethylene, 
and for example in a Monte Carlo study with a lead-hydrocarbon shield [7] the 
authors remark that this "provides effective shielding only for neutrons below 
10 MeV" . 

Thus there is a possibly significant question of "fast" or "higher energy " 
neutrons. It may therefore be useful if we present some simple remarks about 
the effects to be anticipated from them in dark matter detectors and perhaps 
other low background setups. We can do this in a simple and general way terms 
using the 'black disc' model which, for fast neutrons, approximately applies for 
all nuclei. While this can only be regarded as a semi-quantative approximation, 
it offers a simple and general overwiew of the situation. Unfortunately, even 
such an rough discussion is not possible for slow neutrons. 

We should stress that we use the terms "fast" or "high energy " in the sense 
of neutron physics, that is to mean neutrons of a few MeV and above, and of 
course not in the sense of high energy particle physics. 

2 "High energy " neutrons and the "Black disc" 

Above a few MeV neutron-nucleus scattering begins to resemble the "black disc" 
limit. While at lower energies the behavior of the cross section is dominated 
by resonances and shows a complicated structure which varies from nucleus to 
nucleus, at the higher energies a smooth behavior with respect to neutron energy 
sets in, one which can approximately be described by the 'black disc'model. This 
model in its simplest form only involves one parameter, namely R, the radius of 
the nucleus. This permits a simple description, applicable to all target nuclei, 
adequate for our purposes of semi-quantative background estimates. 

In Fig[l] Fig[21 and Fig[3] [5] we show as examples the total experimental 
cross section for the light, medium, and heavy weight nuclei oxygen, calcium, 
and tungsten. One notes that just above the resonance region the total cross 
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section becomes smooth and tends to be above the black disc value, (indicated 
by the grey line using R as given by Eq[3]) and then slowly descend through it 
as the energy increases. 

In the elastic scattering there is a forward diffraction peak at low momentum 
transfer A. This elastic scattering is concentrated in a region characterized by 



where R is the radius of the nucleus. At the same time there is an inelastic 
cross section of about the same magnitude. In the perfect "black disc" limit the 
elastic and inelastic scattering cross sections are equal and of magnitude wR 2 . 
As we shall discuss below, this inelastic scattering could in principle be used as 
a signal for the presence of the fast neutrons. 

2.1 Validity of the simple 'Black disc' model 

The principle features of the simple 'Black disc' model result from the assump- 
tion that the nucleus is totally absorbing for neutrons up to impact parameter 
R, beyond which there is no scattering or absorbtion. This simple assumption 
has the following consequences: 1) the total cross section is made up of equal 
elastic and inelastic parts, adding up to a to t = 2 x irR 2 , 2) these features are 
energy independent 3) the form of the elastic scattering versus angle or momen- 
tum transfer is given by a simple Bessel function expression, (see Eq[7]) where 
in particular the form versus momentum transfer is independent of the neutron 
energy. 

As one sees in the examples of Fig[T] Fig[2] and Fig[3] the energy where the 
smooth 'optical model' behavior sets in varies with the target nucleus. For 
oxygen the energy must be above about 10 MeV, for calcium above 4 or 5 MeV, 
while for tungsten 1 MeV suffices. One also notes from the plots that with our 
value for R the total cross section is somewhat underestimated, particularly at 
the lower energies . 

There is a large literature concerning more sophisticated treatments of the 
optical model [5]. There the simple 'black' assumption is replaced by a partial 
and energy dependent transparency, the 'hard edge' is softened, modifying the 
simple radius assumption, and resonance effects and variation with nucleus may 
be included. While all these points are interesting and important they will not 
significantly alter our general conclusions, particularly our main point, as to be 
discussed in section [3] that the recoil energies implied by the optical model and 
characterized by our parameter E R are in the region anticipated for dark matter 
searches. 

3 Recoil energies 

It is interesting that if we put Eq[T] in the formula for the recoil energy Er 



Ai?~ 1, 



(1) 



A 2 1 
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Er 
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Element 


A 


R(fcrmi) 


E° R (keV) 





16 


3.5 


98 


F 


19 


3.7 


74 


Na 


23 


4.0 


54 


Si 


28 


4.2 


39 


Ar 


40 


4.8 


21 


Ca 


40 


4.8 


21 


Ge 


74 


5.9 


7.7 


I 


127 


7.0 


3.1 


Xc 


132 


7.1 


2.9 


W 


184 


7.9 


1.7 



Table 1: The quantities R, E° R for various nuclei with mass number A. 



one gets recoils in the region of interest for many dark matter searches. M is 
the mass of the nucleus, M w A x 1 GeV and for R we use 

R = A 1 / 3 x 1.4/. (3) 

Calling this typical recoil energy E R one obtains (h, c = 1 units) 

E -=R^M A =^ MeV - (4) 

This is the typical energy characterizing the elastic scattering of the "high energy 
neutron" . In Table 1 we show the value of E° R for some nuclei. The numbers 
run through the values often considered in dark matter searches. 

The differential cross section is a function of the dimensionless parameter 
x = AR, and with these definitions we can write 

(ARf = ^. (5) 

R 

The introduction of the dimensionless variable x, where x 2 = -§£ allows us to 

it 

describe, within the limitations of the model, all nuclei in a universal way in 
terms of their E° R . 

A second interesting point is that in the model this typical recoil energy for 
the "high energy " neutrons depends on the nucleus but not on the energy of the 
neutron. This implies that the whole of the "high energy " neutron spectrum 
contributes in approximately the same way to the recoils. This statement is 
only strictly true insofar as the black disc model is exactly valid, but it will be 
seen that the major features of the model hold approximately over a wide range 
of energy. 

Finally we note the relation between scattering angle and recoil energy 

t-1 ^ (-^neutron) /-, n \ j-^ m neutron/^ n \ 

Er = 2M~ A = M A (1 " C0S9) = 2 ^-*™„^— (1 - cos6) , (6) 
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with Ma the mass of the nucleus. This simple expression will also be true, 
replacing "neutron" with "WIMP" , for the elastic WIMP scattering when the 
mass of the WIMP is small compared to that of the nucleus. 

In the 'black disc' limit the energy spectrum of the recoils is quite simple. 
The diffraction peak leads to an Er spectrum falling off with energy, but at 
least for the lighter nuclei not so steeply as from WIMP scattering. The elastic 
scattering is given by a Bessel function expression [TU], ~ ( Ji< ar^ ) an< ^ 
since Er — A 2 /2Ma, one expects 



dE R dA 2 \ AR J 

Now since (AR) 2 = Er/E r we can introduce the convenient variable x, where 
x 2 = Er/E r and write 



da da ( 2J\{x) 
dER dx 2 



„2 



X 



Er/E° r , (8) 



We include the '2' in the Bessel function expression so it is normalized to 
one at x = 0. For small x, x < 1, the expression then behaves as 



{2Ji(x)Y _« , 2j . , iEr 



V x J 4 "' 4 E° R 



(9) 



For small w this implies a relatively flat energy spectrum. Thus for oxygen, 

between Er = 10 and 40 keV, goes down only by about 8%. Note again 
that it is not necessary to know the neutron energy spectrum to reach this 
conclusion, as long as we have "high energy " neutrons. This is of course quite 
different from the rapid fall-off of WIMP scattering, even on light nuclei. The 
difference originates in the fact that while for fast neutrons the recoil spectrum 
is determined by the radius of the nucleus, for WIMPs it is largely determined 
by the velocity spectrum of the incoming particles, which falls steeply above a 
certain velocity. 

Hence it appears that a characteristic sign of the fast neutron background 
would be a rather flat recoil energy spectrum for the light elements. 

For -if- larger than one the black disc elastic cross section falls rapidly. For 
convenience we show a few values in Table 2. One sees that most of the elastic 

cross section is contained within x 2 = -§#■ < a few 

■ i-i «... 
In Fig|4] we show some data on the angular distribution for oxygen. The 

optical model shape is quite clear, with a well- defined low momentum transfer 

peak and diffractive oscillations, although at this energy the effective value for 

R is slightly smaller than the one we are using. To compare with our variable A 

we have added an axis showing the momentum transfer for 20 MeV neutrons. 
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X 


E R / E° r 


(2j 1 (x)/ x y 








i 


1.0 


1.0 


0.77 


1.6 


2.6 


0.51 


2.0 


4.0 


0.33 


2.6 


7.0 


0.13 


3.0 


9.0 


0.051 


3.6 


13 


0.0028 



Table 2: Some values for the Bessel function expression for the elastic cross 
section in the 'black disc' model, normalized to 1 at zero recoil energy. 

4 Patterns of A Behavior 

Given a certain neutron flux in a detector with different nuclei, the rate will 
vary in a characteristic way with increasing mass number A. In the black disc 
limit the total elastic cross section, that is, integrated over all En, is ttR 2 , and 
so one has 

a ~ R 2 ~ A 2/3 all E R . (10) 

However, most detectors will have a lower energy threshold for detection and 
so will miss part of the diffraction peak. So it is perhaps also useful to consider 
the behavior of -j^- at some fixed Er, where one has 

da _ , MRi / 2^ 2 ^ A7/3 (M^S 2 s= ,E*\ * 



dE R 2 \ x J \ x ) \E° R 

This formula follows from EqJT] and the use of the optical theorem point for 
forward scattering with a purely imaginary amplitude, as one has for the 'black 
disc': 



da 



dA 2 



A=0 167r 



~total ■ (12) 



In terms of the x variable alone, these relations can also be given the simple 
form 

At fixed En, the increase of ~ A 7 / 3 in Eq[TT]will by slowed by the decreasing 
E R , which causes x to move out in the Bessel function expression. Since in the 
'black disc' limit the recoil spectrum is independent of the incident neutron 
energy and the incoming flux is the same for all nuclei in a given experimental 
setup, we may take the rate as approximately proportional to Eq llll In Table 
3 we show, for E R = 20 and 30 keV, how EqlTll and thus the rate varies with A 
for some elements. 
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Element 


A 


neutron 


neutron 


WIMP 


WIMP 






E R =20kcV 


E R =30keV 


M=10 GeV 


M=50 GeV 





16 


1 


1 


1 


1 


F 


19 


1.5 


1.5 


1.3 


1.8 


Na 


23 


2.2 


2.2 


1.6 


3.3 


Si 


28 


3.4 


3.3 


1.8 


6.7 


Ar 


40 


7.0 


6.4 


1.1 


19 


Ca 


40 


7.0 


6.4 


1.1 


19 


Ge 


74 


19 


13 


-0 


93 


I 


127 


20 


5.1 


-0 


200 


Xe 


132 


18 


3.9 


-0 


240 


W 


184 


2.6 


1.6 


-0 


230 



Table 3: Variation of the differential scattering rate per unit energy over various 
nuclei in the 'black disc' limit, Eq[TTJ at E R =20keV and E R =30keV. For 
comparison the same rate for a coherently scattering WIMP at E R — 20 keV 
for masses 10 and 50 GeV is also shown. One notes different patterns of A 
behavior for neutrons and WIMPs. All values are per nucleus and normalized 
to that for oxygen. 

The interplay of the factors in Eq llll leads to a maximum of the rate at fixed 
E R towards the middle of the periodic table. An interesting point is that while 
for lower energy neutrons one expects, on essentially kinematic grounds, the 
scattering to be predominantly on light nuclei such as oxygen, for "high energy" 
neutrons , where the incident energy is high relative to the recoil energy, this is 
no longer true. 

All values in the Table are normalized to that on oxygen. The entries refer 
to the rates per nucleus, so for a given detector account should taken of the 
relative abundance of the nucleus in the target material. 

For comparison we also show the rates for 10 and 50 GeV WIMPs, with 
coherent ~ A 2 interactions, for E R = 20 keV. For WIMPs not having the simple 
coherent interactions, it is not possible to make a simple statement as to the A 
behavior since the variation from nucleus to nucleus will depend on the quantum 
number constitution of the WIMP. Indeed it might be possible to unravel these 
quantum numbers by studying the A behavior [13] . 

E q|13l may be used to estimate the event rate in an acceptance window for a 
given setup and an assumed neutron flux. Bessel function relations can be used 
to estimate what fraction of the elastic cross section irR 2 is in a certain interval. 
Thus, for example, assume a CaW04 detector accepting events on the oxygen 
between lOkeV and 40keV. From Table 1 this correspopnds to x 2 = E R /E° be- 
tween 0.11 and 0.42. On the other hand the total elastic cross section according 
to E q[T3l is proportional to J °° (2Ji(x)/x) 2 dx 2 . Now there is the "orthonormal- 
ity relation" [12] / °° (2J 1 (x)/x) 2 dx 2 = 4, where we have normalized as in Eq[9j 
According to Eq[S]one may take {2J 1 (x)/x) 2 as approximately equal to one in 
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this interval. Thus we have a fraction of the elastic cross section 0.3f/4 = 0.078 
in the interval 10 keV to 40 keV. Since irR 2 — 0.39 b, there is thus a cross section 
of 0.030 b in the interval. 

For a flux of say 1 x 10~ 7 /cm 2 s for the "high energy " neutrons, this leads 
to a rate of 2 x 10~ 3 /kgday for CaWC>4. Some information on underground 
neutron fluxes may be found in Reference P3]. Naturally given the neutron 
flux at a given location, the neutrons must be further be propagated through 
the shielding and materials of the setup in question, and there is the further 
question of neutrons originating in these materials. 

5 Inelastic Processes 

An important difference between neutrons and WIMPs is of course that neu- 
trons are strongly interacting and WIMPs weakly interacting. This implies 
that, depending on the material and geometry of the detector, multiple scatter- 
ing may serve as a signal for neutrons. 

A further significant difference between the MeV neutrons we discuss here 
and the assumed galactic WIMPs is the low energy of the WIMP. The WIMP 
is expected to have a small velocity, typical of objects in the Milky Way, v ~ 
2 x 10~ 3 (c=l units). Then even, say, a mass 50 GeV WIMP only has an energy 
\mv 2 ~ 100/ceF, and a 'light WIMP' with M=10 GeV only - 20fceV. Such 
energies are generally not enough to induce significant inelastic processes on 
most nuclei. 

On the other hand, neutrons with MeVs and above typically have many in- 
elastic reations and the inelastic reactions are expected at close to the same rate 
as the elastic scattering. In the 'black disc' limit' one has (Jeiastic ~ <?inelasUc ~ 
\o 'total ■ Fig[5] j8] for the elastic scattering on oxygen, when compared with Fig[TJ 
shows this is approximately true. 

In experiments with two-signal readout designed to distinguish electromag- 
netic backgrounds from nuclear events, such as CRESST, Edelweiss, CDMS, 
Xenon.. . these inelastic events will typically appear as a high energy signal in 
the nuclear event band. Depending on the nuclei in question, the most com- 
mon inelastic channels will be nuclear excitation, leading to a 7 ray, or nuclear 
breakup. In breakup common nuclear fragments are further neutrons and a 
particles, while for nuclear excitation followed by 7 rays electromagnetic energy 
is present. Therefore inelastic events will generally have a relatively high energy 
deposit, compared to the elastic recoils either from WIMPs or neutrons. De- 
pending on the nuclei present in the detector, these inelastic events could also 
appear on nuclei other than those exhibiting the small recoils. Thus associated 
with a fast neutron flux there should be a characteristic set of events with quite 
different properties from the small elastic scattering recoils. But these should 
be absent for WIMPs. 
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6 Conclusions 



We have pointed out that the difficult-to-shield fast neutrons will give elastic 
nuclear recoils of the same type as sought for in dark matter searches. We 
explain, however, that it should be possible to use evidence internal to the de- 
tector in establishing the presence of a neutron background and so not to be 
only dependent upon standard assumptions and Monte Carlo estimates. We 
have suggested some methods of doing this, including the pattern of recoil en- 
ergies , the rates on different nuclei through the periodic table (see section 4) 
and the presence of inelastic nuclear processes at about the same level as elastic 
nuclear recoils (see section 5). In this respect it is thus desirable to have detector 
systems able to distinguish more than one type of recoiling nucleus |15j . 

For light to medium weight nuclei the neutron induced events should show a 
rather flat spectrum with respect to recoil energy Er. At a fixed recoil energy, 
the variation of the rate with respect to changing the nucleus is different for 
neutrons and WIMPs as explained in connection with Table 3. 

Neutrons can of course have multiple interactions, and fast neutrons as 
opposed to WIMPs can induce inelastic reactions, and at a relatively high rate. 

The great difficulty, as with all aspects of dark matter searches, is that at 
present only a handful of events is available and that a detailed investigation 
of these points requires extensive data, not easy to come by in very low rate 
experiments. 
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Figure 1: The total cross section for neutrons on le O, from the data compi- 
lations of the of the neutron data center NNDC at Brookhavcn National Lab- 
oratory. Above the resonance region the cross section behaves smoothly. The 
'black disc' estimate (grey line), using R = A 1 / 3 !.*!/, is 2 x irR 2 = 0.78 b. 
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Figure 2: The total cross section for neutrons on i0 Ca from 1 to 100 MeV, 
from the NNDC. The black disc value with R = A 1 ' 3 1.4/ is 1.4 b 
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Figure 3: The total cross section for neutrons on 184 W from 1 to 100 MeV, 
form the NNDC. The black disc value with R = A 1 ' 3 1.4/ is 4.0 b 
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9cm. 

H 1 1 > 

36 315 840 E(recoil)/keV 

Figure 4: The angular distribution for elastic neutron-oxygen scattering for 
several energies from ref [II] . The added horizontal axis shows, for the 20 MeV 
case, En for scattering angles 10, 30 and 50 degrees. In the 'black disc' limit 
the curves should be the same for all energies, when plotted against Er instead 
of angle. 
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Figure 5: The total cross section for elastic scattering on le O, from the NNDC. 
As compared to the total cross section in Fig[TJ this elastic cross section includes 
only processes where the nucleus undergoes no breakup or excitation. One notes 
that the simple optical model relation a e iastic ~ \vtotai is approximately correct. 
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